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Abstract 

Ohya and Volovich have proposed a new quantum computation model 
with chaotic amplification to solve the SAT problem, which went beyond 
usual quantum algorithm. In this paper, we generalize quantum Turing 
machine, and we show in this general quantum Turing machine (GQTM) 
that we can treat the Ohya- Volovich (OV) SAT algorithm. 

1 Introduction 

The problem whether NP-complete problems can be P problem has been con- 
sidered as one of the most important problems in theory of computational com- 
plexity. Various studies have been done for many years pQ. Ohya and Volovich 
2 , 3 proposed a new quantum algorithm with chaotic amplification process to 
solve the SAT problem, which went beyond usual quantum algorithm. This 
quantum chaos algorithm enabled to solve the SAT problem in a polynomial 
time glEHl]. 

In this paper we generalize quantum Turing machine so that it enables to 
describe non-unitary evolution of states. This study is based on mathemati- 
cal studies of quantum communication channels E] . It is discussed in this 
generalized quantum Turing machine (GQTM) that we can treat the OV SAT 
algorithm. 

In Section 2, we generalize QTM by rewriting usual QTM in terms of channel 
transformation so that it contains both dissipative and unitary dynamics. In 
Section 3, the SAT problem is reviewed and fundamental quantum unitary gates 
are presented. In Section 4, based on the papers 0][7|, we concretely construct 
the fundamental gates needed for computation of the SAT problem. In Section 
5, we rewrite the total process including a measurement process and amplifier 
process with chaotic dynamics by GQTM. 
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2 Generalized Quantum Turing Machine 

Classical Turing machine(TM or CTM) M cl is denned by a triplet (Q, £,<$), 
where E is a finite alphabets with an identified blank symbol Q is a finite set 
of states (with an initial state go and a set of final states g/) and 6 : Q x £ — > 
Q x £ x {—1, 0, 1} is a transition function. Note that {—1, 0, 1} indicates moving 
direction of the tape head of TM. The deterministic TM has a deterministic 
transition function 6 : QxS — * 2^ xS X {-1,0,1}, that is, S is a non-branching 
map, in other words, the range of 8 for each (q, a) 6 Q x £ is unique. A TM M 
is called non-deterministic if it is not deterministic. 

Quantum Turing machine (QTM) was introduced by Deutsch jS] and has 
been extensively studied by Bernstein and Vazirani [5] .In this section, we intro- 
duce a generalized quantum Turing machine (GQTM), which contains QTM as 
a special case. 

The Hilbert space TL of QTM consists from complex functions defined on 
the space of classical configurations. 

Definition 1 Usual Quantum Turing machine M q is defined by a quadruplet 
M q = (Q, £, TL, U) , where TL is a Hilbert space described below in f2.1p and U is 
a unitary operator on the space TL of the special form described below in $2.2\) . 

Let C=QxExZbe the set of all classical configurations of the Turing 
machine M c i, where Z is the set of all integers. It is a countable set and one has 



Since the configuration C G C can be written as C — (q, A, i) one can say that 
the set of functions {| q,A,i >} is a basis in the Hilbert space TL. Here q E Q, 
i £ Z and A is a function A : Z — * ^. We will call this basis the computational 
basis. 

By using the computational basis we now state the conditions to the unitary 
operator U. We denote the set T = { — 1,0,1}. One requires that there is a 

function ^ :QxExQxSxr^C which takes values in the field of computable 

numbers C and such that the following relation is satisfied: 



p.b,a 

Here the sum runs over the states p S Q, the symbols b £ £ and the elements 
g G r. Actually this is a finite sum. The function A\ : Z — > -ft is defined as 



The restriction to the computable number field C instead of all the complex 
number C is required since otherwise we can not construct or design a quantum 
Turing Machine. 




(2.1) 




(2.2) 
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Note that if, for some integer t £ N = {1, 2, ...} , the quantum state U* |go, A, 0) 
is a final quantum state, i.e. \\EQ(qp)U s \qo,A,0)\\ = 1 and for any s < t, s G N 
one has \\EQ(qp)U s \qo,A, 0)|| = 0, then one says that the quantum Turing 
machine halts with running time t on input A. 

Now we define the generalized quantum Turing machine (GQTM) by using 
of a channel A (see below) instead of a unitary operator U. 

Definition 2 Generalized Quantum Turing machine M gq (GQTM) is defined 
by a quadruplet M gq = (Q, £,7i, A) , where Q and £ are two alphabets, H is a 
Hilbert space and A is a channel on the space of states on H of the special form 
described below. 

Let us explain GQTM in more detailed. GQTM M gq is defined by quadru- 
plet ((5,E,7i,A), where Q is a processor configuration, £ is a set of alphabet 
including a blank symbol and A is a quantum transition function sending a quan- 
tum state to a quantum state. Q and £ are represented by a density operator 
on Hilbert space Hq and Hz, which are spanned by canonical basis {\q} ; q £ Q} 
and {\a) ; a £ £} , respectively. A tape configuration A is a sequence of elements 
of £ represented by a density operator on Hilbert space Hs spanned by a canon- 
ical basis {\A) ; A £ £*} , where £* is the set of sequences of alphabets in £. 
A position of tape head is represented by a density operator on Hilbert space 
Hz spanned by a canonical basis {\i) ;i £ Z}. Then a configuration p of GQTM 
M gq is described by a density operator in H = Hq®Ht.®Hz- Let & (H) be the 
set of all density operators in Hilbert space H. A quantum transition function 
A is given by a completely positive (CP) channel 

A : 6 (H) — > 6 (H) . 

For instance, given a configuration p = J2 k Xk \ipk) (V'fcl ; where ^ Afe = 1, Afe > 
and ipk = |<7fc) ® \Ak) <8> \ik) {qk £ Q, A k £ T>* ,i k £ Z) is a vector in a basis of 
H. This configuration changes to a new configuration p' by one step transition 
as p' = A (p) = J2 k f-k \ipk) (i>k\ with £>fc = 1, > 0. 

One requirement on GQTM M gq = (Q, E,7i, A) is the correspondence with 
QTM. If the channel A in GQTM will be a unitary operator U then GQTM 
M g q = (Q,T,,H,A = U ■ U*) reduces to QTM M q = (Q,Y,,H,U). 

Several studies have been done on QTM whose transition function is repre- 
sented by unitary operator, in which various theorems and computational classes 
in QTM were discussed in PHU|. 

Let us explain how to construct a QTM. Let 8 be a function 

(5:QxExQx£x {-1,0, 1} -> C. 
For any q £ Q, a £ S, it holds 

^2 \8(q,a,p,b,d)\ = 1. 

p£Q,6eS,d£{-l,0,l} 

For any q £ Q, a £ S, q' <?) G Q, a' a) G E, it holds 
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6{q',a',p,b,d)*S(q,a,p,b,d) = 0. 

p£Q,6es,de{-i,o,i} 

Given QTM M q and its configuration p = \ip) ((f \ with \ip) — \q, A, i), after 
one step, this configuration is changed by a transition function 5 as 

A s (\q,A,i)(i,A,q\)= ^ S(q, A(i),p,b,a)S*(q,A(i),p',b',a) 

p.b,(T,p' ,b f ,<T f 

\p, A\,i + a) {i + a ',A h - ,p'\ 

Remark 3 For any q,p S Q, a, b E S, d e { — 1, 0, 1}, /ei <5 (q, a,p, b, d) — {0, 1}, 
then QTM is a reversal TM. 

A transition of GQTM is regarded as a transition of amplitude of each con- 
figuration vector. We categorize GQTMs by a property of CP channel A as 
below. 

Definition 4 A GQTM M gq is called unitary QTM (UQTM, i.e., usual QTM), 
if all of quantum transition function A in M are unitary CP channel. 

For all configuration p = J^n^nPn (S„A„ = 1,A„ > 0), a GQTM M gq is 
called LQTM Mi q if A is affine ; A (J2 n ^npn) — J2 n (Pn) ■ Since a mea- 
surement defined by AmP = ^PkpPk with a PVM {Pk} on Ti is a linear CP 

k 

channel, LQTM may include a measurement process. 

For a more general channel the state change is expressed as 

A(\q,A(i),i){q,A(i),i\)= S(q, A(i),p,b, a,p' ,b' , a') 

p,b,cr,p f ,b' ,a f 

\p, A*, i + a) {p,Af,i + a\ 

with some function 8(q, A(i),p, b, o,p' , b' , a') such that the RHS of this relation 
is a state. 

Thus we define two more classes of GQTM for non-unitary CP channels. 

Definition 5 A GQTM M gq is called a linear QTM(LQTM) if its quantum 
transition function A is a linear quantum channel. 

Unitary operator is linear, hence UQTM is a sub-class of LQTM. moreover, 
classical TM is a special class of LQTM. 

Definition 6 A GQTM M gq is called non-linear QTM (NLQTM) if its quan- 
tum transition function A contains non-linear CP channel. 

A chaos amplifier used in [2 E] is a non- linear CP channel, the details of this 
channel and its application to the SAT problem will be discussed in the sequel. 
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2.1 Computational class for GQTM 

Let us state some language classes which classical Turing machine recognizes. 

Definition 7 The class of languages is in P if its language is recognized by a 
deterministic Turing machine in polynomial time of input size. 

Definition 8 The class of languages is in NP if there is a deterministic Turing 
machine, called the verifier, which recognize languages with some informations 
in polynomial time of input size. Besides, if a language L\ €NP and L\ reduces 
to L12 GNP in polynomial time, a language L\ is NP-complete. 

Definition 9 If languages are accepted by non- deterministic Turing machine in 
polynomial time of input size with a certain probability, this class of languages 
are called the class of bounded probability polynomial time(BPP). 

A NP-complete language is the most difficult one in NP. If there is a poly- 
nomial time algorithm to solve it in the above sense, it implies P=NP. The 
existence of such a algorithm is demonstrated in [21 E] in an extended quantum 
domain, as is reviewed in the next section. In this paper we will show that this 
OV algorithm can be written by GQTM in the sequel section. 

Given a GQTM M gq — (Q, S, S) and an input configuration po — \vin) (vin\t 
i\vin) = ko) ® \T) <8 |0)), a computation process is described as the following 
product of several different types of channels 

Ai o •• - o A t (p ) = Pf = \v f ) (v f \ 

where Ai, • • • , A t are CP channels. Applying the CP channels to an initial state, 
we obtain a final state pf and we measure this state by a projection (or PVM) 

Pf = \Qf) («/|®/e®/z, 

where Is,Iz are identity operators on TLy.,'Hz, respectively. Let p > be a 
halting probability such that 

trn^Hz (PfPf) = P kf) (qf\ ■ 

Then, we define the acceptance (rejection) of GQTM and some classes of 
languages. 

Definition 10 Given GQTM M gq and a language L, if there exists N steps 
when we obtain the configuration of acceptance (or rejection)by the probability 
p, we say that the GQTM M gq accepts (or rejects)L by the probability p, and 
its computational complexity is t. 

Definition 11 A language L is bounded guantum probability polynomial time 
GQTM(BGQPP) if there is a polynomial time GQTM M gq which accepts L with 
probability p > j. 
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Similarly, we can define the class of languages BUQPP(= BQPP), BLQPP, 
BNLQPP(=BGQPP) corresponding to UQTM, LQTM and NLQTM, respec- 
tively. 

In Section 2, it is pointed out that LQTM includes classical TM, which it may 
imply: BPPCBLQPPLCBNLQPPCBGQPP. Moreover, if NLQTM accepts the 
SAT OV algorithm in polynomial time with probability p > i, then we may 
have the inclusion 

NP C BGQPP 

We will discuss this inclusion in Sec. 4 by constructing GQTM which accepts 
the SAT OV algorithm. 

3 SAT Problem 

Let X = {xi, . . . , x n } , n G N be a set. a; fc and its negation x/. (k = 1, . . . , n) 
are called literals Let X = {xi, . . . ,3^"} be a set, then the set of all literals is 
denoted by X' = X U X = {x\ , . . . , x n ,xi, . . . , x^}. The set of all subsets of X' 
is denoted by T {X') and an element C G J-(X') is called a clause. We take a 
truth assignment to all variables Xk- If we can assign the truth value to at least 
one element of C, then C is called satisfiable. When C is satisfiable, the truth 
value t (C) of C is regarded as true, otherwise, that of C is false. Take the truth 
values as "true <->l, false <->0". Then Cis satisfiable iff t (C) = 1. 

Let L = {0, 1} be a Boolean lattice with usual join V and meet A, and t (x) 
be the truth value of a literal x in X. Then the truth value of a clause C is 
written as t (C) = V xe ct (x). 

Moreover the set C of all clauses Cj (j = 1, 2, • • • , m) is called satisfiable iff 
the meet of all truth values of Cj is 1; t(C) = Af =l t{Cj) = 1. Thus the SAT 
problem is written as follows: 

Definition 12 SAT Problem: Given a Boolean set X = {x\, • • • , x n }and a set 
C = {Ci, • • • , C m } of clauses, determine whether C is satisfiable or not. 

That is, this problem is to ask whether there exists a truth assignment to 
make C satisfiable. It is known in usual algorithm that it is polynomial time to 
check the satisfiability only when a specific truth assignment is given, but we 
can not determine the satisfiability in polynomial time when an assignment is 
not specified. 

In 0] we discussed the quantum algorithm of the SAT problem, which was 
rewritten in [7] with showing that the OM SAT-algorithm is combinatorial. In 
|2 13] it is shown that the chaotic quantum algorithm can solve the SAT problem 
in polynomial time. 

Ohya and Masuda pointed out 0] that the SAT problem, hence all other 
NP problems, can be solved in polynomial time by quantum computer if the 
superposition of two orthogonal vectors |0) and |1) is physically detected. How- 
ever this detection is considered not to be possible in the present technology. 
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The problem to be overcome is how to distinguish the pure vector |0) from the 
superposed one a |0) + /3 |1) , obtained by the OM SAT-quantum algorithm, if 
f3 is not zero but very small. If such a distinction is possible, then we can solve 
the NPC problem in the polynomial time. In El it is shown that it can be 
possible by combining nonlinear chaos amplifier with the quantum algorithm, 
which implies the existence of a mathematical algorithm solving NP=P. The al- 
gorithm of Ohya and Volovich is not known to be in the framework of quantum 
Turing algorithm or not. This aspect is studied in this paper. 



3.1 Quantum computation 

In this subsection, we review fundamentals of quantum computation (see, for 
instance, JD). Let C be the set of all complex numbers, and |0) and |1) be the 
two unit vectors (q) and (^) , respectively. Then, for any two complex numbers 
a and (3 satisfying |a| 2 + |/3| 2 = 1, a |0) +(3 11) is called a qubit. For any positive 
integer N, let H be the tensor product Hilbert space defined as (C 2 ) and let 
■T | ei) ; < i < 2 7V_1 } be the basis whose elements are denoted as 



eo) 


= |0)<5 


§> |o> - - - 


<8> |0) EE 


|o,o, • - 


• ,0), 


ex) 


= |1>« 


5|0>- 


® |0) EE 


|i, o, - - 


• ,0), 


e 2 > 


= |0>« 


5|1>- 


® |0) EE 


|0, 1, - • 


• ,0), 



|e 2W _ 1 > = |l)®|l)-..®|l) = |l,l J .-.,l>. 

For any two qubits \x) and \y), \x,y) and \x N ) is defined as \x) <g \y) and 
\x) g> ■ ■ ■ g> \x) , respectively. 

v 

N times 

The usual (unitary) quantum computation can be formulated mathemati- 
cally as the multiplication by unitary operators. Let Unot,Ucn and Uccn be 
the three unitary operators defined as 

U N ot = |1) (0| + |0) (1| , 

U CN = |0) (0|®J+|1)(1|®I^ot, 

Uccn = |0) (0| ® I® / + |1) (1| g> |0) (0| ® 1+ |1) (1| g> |1) (1| g> U NOT - 

Unot,Ucn and Uccn represent the NOT-gate, the Controlled-NOT gate and 
the Controlled-Controlled-NOT gate, respectively. Moreover, Hadamard trans- 
formation H is defined as the transformation on C 2 such as 

ff|o) = -^(|o) + |i)), tf|i) = -L(| )-|i». 

These four operators Unot, Ucn, Uccn and H are called the elementary gates 
here. For any k G N, (k) denotes the fc-tuple Hadamard transformation 

on (C 2 )^ defined 
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2 k ~ 1 

U H ] (*) K) = ^2 d°) + l 1 ))^ I """) = Wfi E ^ ® 1°^) ■ 

The above unitary operators can be extended to the unitary operators on 

U ( N N J T (u) ee I^- 1 ® (|0) <1| + |1) (0|) J«"-«-i 
C/W (u, «) ee J®"" 1 <g> |0) (0| ® J®^""" 1 + g, |l) (i| 

I/^jv («> v ' w ) = ® 1°) (°l ® i"®^"" 1 + I®"" 1 ® |1> (1| 

® g | ) (o| $ 

+ I® 11-1 ® |1) (1| <g> j®"-"- 1 8 |i) (i| 8 

where u, i> and w be a positive integers satisfying l<u<v<w<N. 

Furthermore we have the following three unitary operators Uand, Uor and 
Ucopy j called the logical gates; (see [7]) 

Uand = £2, £i A £ 2 ) (ei, £ 2 , 0| + |£i, £ 2 , 1 - £i A £ 2 ) (£i, £ 2 , 1|} 

ei,e 2 e{0,l} 

= |0,0,0) (0,0, 0| + |0,0,1) (0,0, 1| + |1,0,0) (1,0, 0| + |1,0,1) (1,0, 1| 
+ |0,1,0) (0,1,0| + |0,1,1) (0,1,1| + |1,1,1)(1,1,0| + |1,1,0)(1,1,1|. 



U OR = ^2 {|£i,£ 2 ,£i V e 2 ) (ei,£2,0| + |ei,e 2 , 1 - £1 V e 2 ) (ei,£2, 1|} 

El, E26{0,1} 

= |0,0,0) (0,0, 0| + |0,0,1) (0,0, 1| + |1,0,1) (1,0,0| + |1,0,0) (1,0, 1| 
+ |0, 1, 1) (0, 1, 0| + |0, 1, 0) (0, 1, 1| + |1, 1, 1) (1, 1, 0| + |1, 1, 0) (1, 1, 1| . 



Ucopy = {l £ i> £ i) ( £ i>°l + l £ i> 1 ~ £ i) ( £ i> !|} 

£lG{0,l} 

= |0,0) (0,0| + |0, 1) (0, 1| + |1, 1) (1,0| + |1,0) (1, 1| . 

We call Uand, Uor and Ucopy , AND gate, OR gate and COPY gate, respec- 
tively, whose extensions on (C 2 )® N are denoted by U a ^ d ,Uq2 and Uqq PY , 
which are expressed as 
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U { A N J D (u,v,w)= ^ 8ti_1 ® |ei> <eil Z®"-"- 1 (g> |e a > <e 2 | 

ei,e 2 £{0,l} 

jVw-v-u-l g | £l A £2 ) (J | 7 ®jV- w -„-„ + 

J®"" 1 |e x ) (eil/^-"- 1 ^ 

|e 2 ) (e 2 | ® |l - ei A e 2 ) (1| jSW-f-»-« 

U^{u,v,w)= ]T ^ U " 1 ®ki)(£i|^""" 1 ®N)<£ 2 | 

ei,e a e{0,l} 

j®iu-t>— u-1 g | £i v £ ^ J0 W- w _„-„ + 

J® u - 1 ®|ei)(ei|r®"- u - 1 ® |e 2 ) (e a | 

4oW («.«) = E 7 ®" _1 ( £ il J®"""- 1 ki> (0| 
eie{o,i} 

+ I®"" 1 |ei> (ei| /»«-«- 1 |i - ei ) (l| j»JV-«-». 

where u, u and u> are positive integers satisfying l<u<v<w<N. These 
operators can be written, in terms of elementary gates, as 

U OR ( u > v > w ) = U CN ( u > w ) ■ U CN ( w : w ) ■ U CCN ( u > v > w ) > 
U AND ( U > w > w ) = ( U ' W ' W ) > 

4 SAT Algorithm 

In this section, we explain the algorithm of the SAT problem which has been 
introduced by Ohya-Masuda |1| and developed by Accardi-Sabbadini [7|. The 
computation of the truth value can be done by by a combination of the unitary 
operators on a Hilbert space Tl, so that the computation is described by the 
unitary quantum algorithm. The detail of this section is given in the papers 
011313 ?]i so we will discuss just the essence of the OM algorithm. Throughout 
this section, let n be the total number of Boolean variables used in the SAT 
problem. Let C be a set of clauses whose cardinality is equal to m. Let TL = 
(C 2 ) n+IJ " be a Hilbert space and |«o) be the initial state \vq) — |0", M , 0) , 
where fi is the number of dust qubits which is determined by the following 
proposition. Let £/g be a unitary operator for the computation of the SAT: 

4 n) \v ) = ^= £ \ei,x»,t ei (C))=\v f ) 

v 2 i=0 
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where x^ L denotes the /i strings in the dust bits and t ei (C) is the truth value of 
C with ei . In HE), U^ l) was constructed. 

Let {sk\ k = 1, . . . , m} be the sequence defined as 

s\=n + l, 

s 2 = Si + card (C x ) + 6 ltCard ( Cl ) - h 

Si = Sj_i + card (d-i) + 6 1 ^ card ( C ._ l) , 3 < i < m, 

where card{Ci) means the cardinality of a clause Cj. And let define s/ as 

Sf = s m - 1 + card(C m ) + S liCard (c m )- 

Note that the number m of the clause is at most 2n. Then we have the following 
proposition and theorem [?] . 

Proposition 13 For m > 2, the total number of dust qubits \i is 
/j, = Sf — 1 — n 

m 

= ^ Card (Cfc) + &l,card(C k ) - 2 - 
fc=l 

Determining fi and the work spaces for computing t (Cfc), we can construct 
Uq concretely. We use the following unitary gates for this concrete expression: 



7t( x ) (h\ - J U and ( s fe+i _ Mfe+2 - 2,s fe+2 - 1) , l<fc<m-2 

u AND \ K ) - \ jAx) , , in; i 



Cob ( fc ) = (^4, s fc - card (Cfc) - 1, s fe - card (C fc ) - 2) flgg (Z 3 , Sfc, s fc + 1) U%1 {h,h, s k ) , 

f UqI (u, v, w) , i„e C fe 

(u, to) = I U { *l T (u) • u£g (u, w, u;) • [/^ T («) , x„€C t 

> U NOT ( u ) ■ U NOT ( v ) ■ U OR ( M ' U ' W ) ' U NOT ( u ) ' U NOT ( v ) J X U ,X V £ C k 

where ii, ?2 5 Z3, ^4 are positive integers such that ie z G Cfc or x z G Cfc, (z = ii, . . . , I4). 

Theorem 14 The unitary operator U^ n \ is represented as 

[/(») = C /(^+ 1 ) ( m _ l) . f/^ 1 ) ( m _ 2 ) L/<"+£ +1) (1) 

' u or^ +1) M ■ U ( o^ +1) (m - 1) C/ " fl +A1+1) (1) ■ C/^ +fl+1) (n) . 
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4.1 The resulting state in the SAT algorithm 

Applying the above unitary operator to the initial state, we obtain the final state 
p. The result of the computation is registered as \t (C)) in the last section of the 
final vector, which will be taken out by a projection P n +u,,i = I® n+ ^ <g) |1) (1| 
onto the subspace of H spanned by the vectors \e n ,e^, 1).. 
The following theorem is easily seen. 

Theorem 15 C is SAT if and only if 

Pn+„M n) \V ) ? 

According to the standard theory of quantum measurement, after a mea- 
surement of the event P n +ij,,i, the state p = \vf >< v/\ becomes 

Pn+u,lpPn+u,l _ 

P — 5 = : P 

1 rpP n +n,\ 

Thus the solvability of the SAT problem is reduced to check that p' ^ 0. The 
difficulty is that the probability 



2 _ \T(Cq) 



TrpP n+ ^i = \\P n +v>,i \vf) || = 

is very small in some cases, where |T(Co)| is the cardinality of the set T(Co), of 
all the truth functions t such that t(Co) = 1. 

We put q = •i/Jr with r = \T(Cq)\ . Then if r is suitably large to detect 
it, then the SAT problem is solved in polynomial time. However, for small r, 
the probability is very small so that we in fact do not get an information about 
the existence of the solution of the equation t(Co) = 1, hence in such a case we 
need further deliberation. 

Let go back to the SAT algorithm. After the quantum computation, the 
quantum computer will be in the state 



\v f ) = ^Jl-q 2 |<A>) ® |0) + q \<pi) <8 |1) 



where \<pij and \<pa) are normalized n (=n + /i) qubit states and q = ^Jr/2 n . 
Effectively our problem is reduced to the following 1 qubit problem: The above 
state \vf) is reduced to the state 



|V;) = v / T^|0) + g|l), 

and we want to distinguish between the cases q — and q > (small positive 
number). Let us denote the correspondence from po = l^o) (^o| with p by a 
channel A/; p = A/po- 

It is argued in that quantum computer can speed up NP problems 
quadratically but not exponentially. The no-go theorem states that if the inner 
product of two quantum states is close to 1, then the probability that a mea- 
surement distinguishes which one of the two is exponentially small. And one 
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may claim that amplification of this distinguishability is not possible in usual 
quantum algorithm. At this point we emphasized [3] that we do not propose 
to make a measurement which will be overwhelmingly likely to fail. What we 
did is a proposal to use the output of the quantum computer as an input 
for another device which uses chaotic dynamics. The amplification would be not 
possible if we use the standard model of quantum computations with a unitary 
evolution. However the idea of the paper is different. In \Bj it is pro- 

posed to combine quantum computer with a chaotic dynamics amplifier. Such a 
quantum chaos computer is a new model of computations and we demonstrate 
that the amplification is possible in the polynomial time. 

One could object that we do not suggest a practical realization of the new 
model of computations. But at the moment nobody knows of how to make a 
practically useful implementation of the standard model of quantum computing 
ever. It seems to us that the quantum chaos computer considered in [2] deserves 
an investigation and has a potential to be realizable. 

4.2 Chaotic dynamics 

Various aspects of classical and quantum chaos have been the subject of numer- 
ous studies f and ref's therein). Here we will briefly review how chaos can 
play a constructive role in computation (see EI for the details). 

Chaotic behavior in a classical system usually is considered as an exponen- 
tial sensitivity to initial conditions. It is this sensitivity we would like to use 
to distinguish between the cases q = and q > discussed in the previous 
subsection. 

Consider the so called logistic map which is given by the equation 

x n+ i = ax n (l - x n ) = g(x), x n € [0, 1] . 

The properties of the map depend on the parameter a. If we take, for example, 
a = 3.71, then the Lyapunov exponent is positive, the trajectory is very sensitive 
to the initial value and one has the chaotic behavior [jJJ . It is important to notice 
that if the initial value xq = 0, then x n = for all n. 

The state \tp) of the previous subsection is transformed into the density 
matrix of the form 

p = q 2 Pi + (1 - q 2 ) Po 

where Pi and Po are projectors to the state vectors |1) and |0) . One has to 
notice that Pi and Po generate an Abelian algebra which can be considered as 
a classical system. The density matrix p above is interpreted as the initial data, 
and we apply the channel A = Ac a due to the logistic map as 

a (1 + 9 (p) 03) 

A CA (P) = ^ ' 

where / is the identity matrix and 03 is the z-component of Pauli matrices. 

Pk = A&A (P) 
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To find a proper value k we finally measure the value of 173 in the state pk such 
that 

M k = trp k a 3 . 

We obtain [3] 
Theorem 16 

P k = iI + 9 Y )a3) > ™dM k = 9 \c?). 

Thus the question is whether we can find such a A: in polynomial steps of 
n satisfying the inequality M/. > \ for very small but non-zero q 2 . Here we have 
to remark that if one has q — then p — Pq and we obtain M& = for all k. If 
q 7^ 0, the chaotic dynamics leads to the amplification of the small magnitude q 
in such a way that it can be detected. The transition from p to p fe is nonlinear 
and can be considered as a classical evolution because our algebra generated by 
Po and Pi is abelian. The amplification can be done within at most 2n steps due 
to the following propositions. Since g k (q 2 ) is Xk of the logistic map Xk+i — g{xk) 
with xq = q 2 , we use the notation Xk in the logistic map for simplicity. 

Theorem 17 For the logistic map x n+ \ = ax n (1 — x n ) with a S [0,4] and 
x$ E [0, 1], let xq be 7^- and a set J be {0, 1,2, ... ,n, ... , 2n}. If a is 3.71, then 
there exists an integer k in J satisfying x k > \- 

Theorem 18 Let a and n be the same in above theorem. If there exists k in J 
such that Xk > \, then k > log " 3 ~ 7 \_ 1 ■ 

Corollary 19 If xq = ^ with r = \T (C)\ and there exists k in J such that 
Xk > \i then there exists k satisfying the following inequality if C is SAT. 

5 , 

4 (» - !) • 

From these theorems, for all k, it holds 

f = iff C is not SAT 
k \ >0 iff C is SAT 

5 SAT algorithm in GQTM 

In this section, we construct a GQTM for the OV SAT algorithm. The GQTM 
with the chaos amplifier belongs to NLQTM because the chaos amplifier is 
represented by non-linear CP channel. The OV algorithm runs from an initial 
state po = |i>o) (t>o| to ~p k through p = \vf) (vf\ . The computation from po = 
l^o) ( v o\ to p = \vf) (vf \ is due to unitary channel Ac = Uc • Uc, and that 
from p = \vf) (vf \ to py is due to a non-unitary channel A Cj4 o A/, so that all 



1 - log 2 



log 2 3.71 - 1 



< k < 
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computation can be done by Aq A oA;o Ac, which is a completely positive, so 
the whole computation process is deterministic. It is a multi-track (actually 4 
tracks) GQTM that represents this whole computation process. 

A multi-track GQTM has some workspaces for calculation, whose tracks 
are independent each other. This independence means that the TM can operate 
only one track at one step and all tracks do not affect each other. Let us explain 
our computation by a multi-track GQTM. The first track stores the input data 
and the second track stores the value of literals. The third track is used for 
the computation of t (Cj) , (i = 1, • • • , m) described by unitary operators. The 
fourth track is used for the computation of t (C) denoting the result. The work 
of GQTM is represented by the following 8 steps: 

• Step 1 : Store the counter c = in Track 1. Calculate [| (n — 1)1 + 1, we 
take this value as the maximum value of the counter. Then, store it in 
Track 4. 

• Step 2 : Calculate c + 1 and store it in Track 4. 

• Step 3 : Apply the Hadamard transform to Track 2. 

• Step 4 : Calculate t (Ci) , • • • t (C m ) and store them in Track 3. 

• Step 5 : Calculate t (C) by using the value of the third track, and store 
t (C) in Track 4. 

• Step 6 : Empty the first, second and third Tracks. 

• Step 7 : Apply the chaos amplifier to the result state obtained up to the 
step 6. 

• Step 8 : If c = [f (n - 1)] + 1 or GQTM is in the final state, GQTM halts. 
If GQTM is not in the final state, GQTM runs the step 2 to the step 8 
again. 

Let us explain the above steps for unitary computation (OM algorithm; i.e., 
up to the steps 6 above) by an example. Let the number of literals be n and 
that of clauses be m. Then the language is represented by the following strings 

O n XYlZiCsG{Ci)C E , 

where 

G (Ci) = £\£2 ■ ■ ■ £ n Y~£i£~2 ■ ■ ■ £n 

_ J o k i h 

£k -\i keii 
__/ o kii[ 
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and X, Cs, Y, Ce are used as particular symbols of clauses. For example, given 
X = {1,2,3}, C = {d,C 2 ,C 3 }, & ='({1,2}, {3}), C 2 = ({3},{2}),C 3 = 
({1} , {2, 3}) , the input tape will be 



OOOXCsIIOYOOICeCsOOIYOIOCeCsIOOYOUCe 



First, our GQTM applies DFT to a part of literals on the track 2. The 
transition function for DFT is written by the following table. Put the vector in 
Hq by q. instead of \q.) and denote the direction moving the tape head by R 
for the right and L to the left (Note that O is the starting position). 



The tape head moves to the right until it reads a symbol Cs- When the 
tape head reads Cs, GQTM increases a program counter by one, while moves 
to the right until it reads 1. Then GQTM stops increasing the counter and the 
tape head moves to the top of the tape. According to the program counter, the 
tape head moves to the right as reducing the counter by one. When the counter 
becomes zero, GQTM reads the data and calculates OR with the data in the 
track 2, then GQTM writes the result in the track 3. GQTM goes back to the 
top of the track 1 and repeats the above processes until it reads Y. 

When GQTM reads Y, it calculates OR with the negation and repeats the 
processes as above. When it reads Ce, it writes down fc k in the track 3 and 
clean the workspace for the next calculation. Then GQTM reads the blank 
symbol and it begins to calculate AND. The calculation of AND is done on 
the track 4. GQTM calculates them as moving to the left because the position 
of the tape head is at the end of the track 3 when the OR calculation is finished. 
Then the result of the calculation is showed on the top of the track 4. 

The transition function of OR calculation is described, similar as classical 
TM, by the following three tables: 



# 1 

q a ,0,R q a A,R 

q a ,0,R 9a,l,-R 

Qf,#,R -j^Qb,0,L+ -^=ip b ,l,L ^=q b ,0,L- ^=q b ,l,L 



X 



q bl X,L 



15 








l 


C s 


X 


Y 


# 


9o 


q a , 0, R 




qbi Cs, R 








9a 


q a ,0,R 


q a , 1,R 


qb,C s ,R 


q a , X, R 


quR' Y ' R 


qAND, #, N 


Qb 


qbs,0,R 


<7c,i,0,L 










Qb,l 


<?6,2,0,i? 


9c,i,0,£ 










9b, fc 


qb,k+i,o,R 


<Zc,fe+i, 0, L 










Qb,n-1 


qb,n,0,R 


qc,n, 0, L 










9b, n 










qmi , Y, R 




9c, 1 


q c ,i,o,L 


9c,i, 1,L 


9c, 1 1 Cs, L 


q c ,i,x,L 




q<i,i,#,R 



qc,n qc,n,0,L q c n ,l,L q c ^ n ,Cs,L q c n ,X,L 9d,n,#, -R 

9d,i 9t2,o,0,iV gr t2)1 ,l,JV 

9d,2 9d,i,0,-R 9d,i,l,-R 

9d,fc 9d,fe-i,0,-R 9d,fc-i, 



qa.n 


9d,n-l,0, i? 


9d,rt-l,0, i? 











1 


y 


C s 


9oK 


9 s ,i,,0,i? 


9ft,i, 0, i 






9e 


9e,0,i? 






q e ,C s ,R 


99 


9 s ,i,0,i? 


9ft, i,0, L 






9 s ,i 


q 9 ,2,o,R 


9ft,2,0,i 






9g,fc 


9g,fc+l, 0, i? 


9ft,fc+i,0, L 






9ff,n-l 


9g,n, 0, R 


9ft, n, 0, i 






9s, n 










9ft, 1 


9/1,1,0, £ 


9ft,i,l,i 


9/i,i, y, £ 


9ft, 1, Cs, £ 


9ft, n 


9ft. n, 0, ^ 


9ft,n, 1, i 


qn,n,y,L 


qh,n,Cs,L 


9i,i 


9*2,1, 0,JV 


9*2,0,1, 






9»,2 


9i,i,0,-R 


9i,l, 1, R 






9i,fc 


9i,fe-l,0, R 


9i,fe-l, 1, R 






qi,n 


9i,n-l,0, i? 


9i,n-l, 0, -R 






qj 


9j,0,i 




qj,0,L 


qt2,a,0,N 



qj,0,L 

qh,i,X,L 9»,i,#, i2 

qh,n,X,L q iy n,#,R 
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1 # 

9*3,0 5*3,0,0, R q t3tl ,l,R q a ,0,N 

<7t3,i 9t3,i,0, R g t3 ,i,l,i? <7a,l,^V 

qt3,a 9t4,o,0, N q t 4,sA,N 

qt3,b Qt3,b,#,L <lt3,b,#,L 9t3,c,#, R 

qt3,c q a ,0,N 

The transition function of AND calculation is described by the following 
table: 

1 # 

944,0 te,6,0, R 

qtiA 9t3,b,l,-R 
qAND qti,a,#,L 
qt4 : a <Zt4,a> £ 9i4,6,#,i 9t4,c> -R 

quM qt4,a,if,L <7i4,&,#,£ qt4,d,#,R 
Qu,c 3/,0,£ 

Let g6 be the processor state of GQTM after the step 6 and Tj, i = 1, . . . , 4 
be the strings of the «-th track. Then the OM algorithm showed that the com- 
putation of the SAT problem of the example given above gives us the resulting 
state ps expressed as 

pe = q 2 \ga) (q 6 \ ® |r 1 ,T 2 ,r 3 ,r 4 (i)) (^,^^3^4(1)1 <g> \0) (0\ 

+ {l-q 2 ) \q e ) (©I® \T U T 2 ,T3,T 4 (0)) (T u T 2 ,T3,T A (0)\®\O) (0\, 

where T 4 (1) (resp. T 4 (0)) indicates that the value in the track 4 is 1 (resp. 0). 

Next step, as the three tracks (1,2,3) can be empty, we can apply the chaos 
amplifier to the above p§ in the following manner: 

The transition function of the step 7 denoted by the chaos amplifier is for- 
mally written as 

A& (P6) = 9 k {q 2 ) \<tr) (<tr\ ® |T 4 (1)> (T 4 (1)| ® \0) (0\ 

+ (l-g k (q 2 )) \q 7 ) (q 7 \ <8> |T 4 (0)) (T 4 (0)| ® |0) (0| 

where g is the logistic map explained in Section 4.2. According to El GQTM 
halts in at most [| (71 — 1)] steps with the probability p > 5, by which we can 
claim that C is SAT. 



5.1 Computational complexity of the SAT algorithm 

We define the computational complexity of the OV SAT algorithm as the prod- 
uct of T Q (ll^ and T CA (n) ,where T Q (Uq 1 ^ is the complexity of unitary 

computation and Tqa (n) is that of chaos amplification. 

The following theorem is essentially discussed in ^1 0] . 
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Theorem 20 For a set of clauses C and n Boolean variables, the computational 
complexity of the OV SAT algorithm including the chaos amplifier, denoted by 
T (C,n), is obtained as follows. 

Tgqtm (C, n) = T Q (V< n) ) T CA (n) - O {poly (n)) , 

where poly (n) denotes a polynomial of n. 

The computational complexity of quantum computer is determined by the 
total number of logical quantum gates. This inequality implies that the compu- 
tational complexity of SAT algorithm is bounded by O (n) for the size of input 
n while a classical algorithm is bounded by O (2™) . 
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